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SYSTEM AND METHOD FOR NON-DESTRUCTIVE FUNCTIONAL IMAGING AND 
MAPPING OF ELECTRICAL EXCITATION OF BIOLOGICAL TISSUES USING 
ELECTROMAGNETIC FIELD TOMOGRAPHY AND SPECTROSCOPY 

CROSS-REFERENCE TO RELATED APPLICATION 

[0001] This application is entitled to the benefit of, and claims priority to, 

provisional U.S. Patent Application Serial No. 60/429,272 filed November 26, 2002 and entitled 
"SYSTEM AND METHOD FOR NON-DESTRUCTIVE FUNCTIONAL IMAGING AND 
MAPPING OF ELECTRICAL EXCITATION OF BIOLOGICAL TISSUES USING 
ELECTROMAGNETIC FIELD TOMOGRAPHY AND SPECTROSCOPY," the entirety of 
which is hereby incorporated by reference. 

BACKGROUND OF THE PRESENT INVENTION 
Field of the Present Invention 

[0002] The present invention relates generally to electromagnetic field tomography 
and spectroscopy, and in particular, to the non-invasive fxmctional imaging, detection and 
mapping of electrical excitation of a biological tissue with the help of electromagnetic field 
tomography and spectroscopy using a sensitive material (solution), injected into the biological 
tissue or in the circulation system, that is characterized by having dielectrical properties that are a 
function of the electrical field generated by biological excited tissue. The invention includes 
several versions of the systems differentiated on the basis of multiple frequency, polarization and 
type of sensitive material (solution) utilization. Further, the invention includes computer 
implemented software specifically configured and tailored for the system and method for non- 
invasive detection and mapping of electrical excitation of biological tissue with the graphical and 
three-dimensional tomographic imaging interface. 

Background 

[0003] It has long been known that within the electromagnetic spectrum, biological 
tissues have different electrical/dielectrical properties, and consequently, visual images of such 
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tissues may be produced based on these properties. For example, it is known that the dielectrical 
properties of tissues with high (muscle) and low (fat and bone) water content are significantly 
different. During the last decade the changes in the dielectrical properties of tissues caused by 
various physiological and pathological alterations have been intensively studied. It has been 
demonstrated that dielectrical properties of malignant tumors and normal tissues are different in 
breast, lung, colon and liver. It has also been demonstrated that ischemia, infarction and hypoxia 
change the dielectrical properties of myocardial tissue. The amount of those differences (i.e., the 
contrast in their dielectrical properties) varies with frequency, type of tissue and the presence and 
type of disease, and the magnitude of variation may range from a few percentage points up to 4-5 
times between the normal tissue and the diseased tissue. These examples demonstrate a high 
potential for the use of electromagnetic tomography in biomedical applications. 

[0004] As a result, microwave tomography has been in the scope of interest of 
research groups for several years. For example, research by at least some of the present 
inventors has resulted in U.S. Patent Nos. 5,715,819, 6,026,173 and 6,333,087 for microwave 
tomographic and spectroscopic systems and methods for detection of physiological and 
pathological conditions of biological tissues and physiological imaging of such tissues. The 
entirety of each of these patents is incorporated herein by reference. 

[0005] Unfortunately, the research and development of this technology for 
biomedical applications has also met with significant difficulties. One such difficulty is the high 
attenuation of electromagnetic fields within the body. Attenuation is less at lower frequencies, 
but unfortunately, lower frequencies also result in lower spatial resolution. The compromise 
between attenuation and spatial resolution forms a frequency optimum for microwave imaging. 
J.C. Lin theoretically estimated that the frequency spectrum from 2 GHz to 8 GHz is the 
opfimum for microwave imaging of biological tissue. Our estimations suggest that microwave 
imaging of whole scale biological objects with reasonable acquisition time and spatial resolution 
of 6-8mm can be performed at frequencies near 1 GHz. Of course, experimentally achieved 
spatial resolution cannot compete with the spatial resolution achieved in X-ray imaging, simply 
because of the large difference in wavelength. 
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[0006] However, the possibility of imaging physiological and pathological 
conditions of tissues, highlighted earlier, makes this technology promising. For example, it has 
been determined that microwave tomography and spectroscopy are capable for detection of 
changes in myocardial blood supply, tissue hypoxia, myocardial ischemia and infarction, i.e. 
functional imaging. Research by at least some of the present inventors has experimentally 
proved that a tissue's dielectrical properties are a sensitive indicator of its functional and 
pathological conditions and the degree of such changes is large enough to be tomographically 
imaged. 

[0007] Fig. 1 is a graphical representation of the changes of myocardial e" 
following short time 20%, 40%, 60% and 100% blood flow reduction. The changes of 
myocardial e" following short time 20%, 40%, 60% and 100% blood flow reduction. 
Summarized data for group of seven canines are presented as normalized on the baseline values. 
Data for three frequencies (0.2 GHz, 1.1 GHz, and 6.0 GHz) are expressed as mean +/- SD. 

[0008] Figs. 2A & 2B are graphical illustrations of the spectral changes in 
myocardial permittivity e' and resistance p, respectively, during 10% hypoxia. The percent 
difference from the mean baseline data, summarized for the group of 7 canines, is shown. The 
bar graph inserted into the bottom right of Figure 2A represents group averaged changes in 
arterial blood pH and pOa. 

[0009] Figs. 3A and 3B are graphical illustrations of the spectral changes in 
myocardial permittivity s' (A) and resistance p (B) during 2 hours acute ischemia. The percent 
difference from the mean baseline data, summarized for a group of 6 canines is shown. 

[0010] Fig. 4 is a graphical illustration of the changes in myocardial dielectric 
properties (&@) for 2-week-old canine myocardial infarction. Summarized data for a group of 
five canines are presented as mean percent change +/- SD in dielectrical values from normal 
zones of the infarcted hearts. The values are compared with fresh tissue from a 10-year-old 
human post-infarction aneurysm. 
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[0011] Figs. 5 A & 5B are reconstructed electromagnetic tomographic images of 
excised canine heart (longitudinally view through the long axis base to apex for s' - top and 
transversal view through an area with significant infarction injury for £" - bottom) together with 
anatomical slices. The frequency is 1 GHz, and the scales are in centimeters. 

[0012] However, microwave tomography and spectroscopy do not appear to be 
capable of detecting changes in the dielectrical properties of myocardium, caused by a spread of 
electrical excitation in tissue. Preliminary studies, conducted in tissue bath using cardiac excited 
tissue and electro-mechanical uncoupling pharmaceutical agents, have indicated that the 
dielectrical properties of cardiac excited tissue change during the excitation cycle. The exact 
degree of such changes is unknown at present time, but it appears to be less than 1 percent. This 
is a relatively small variation in tissue dielectrical properties to be reliably reconstructed using 
modem electromagnetic tomographic technologies. Further, it is difficult to acquire the 
necessary data during the short period of time available during the circulation cycle. 

[0013] Thus, previous approaches to localizing the origin of such phenomenon as 
cardiac arrhythmias have depended on one of three principal techniques: catheter mapping, 
electrical excitation mapping during cardiac surgery, or body surface mapping of electrical 
potentials and magnetic fields. Each of these techniques has limitations. For example, catheter 
mapping and excitation mapping during surgery are inherently invasive, provide only limited 
access, and are time sensitive. On the other hand, body surface mapping can be performed in a 
non-invasive, low-risk manner, but with such poor definition that the data is generally considered 
xmsuitable for directing subsequent therapy. Thus, a need exists for a non-invasive system by 
which electrical excitation of a biological tissue may be reliably detected and mapped. 

[0014] Another difficulty faced in the use of microwave tomography for the purposes 
described hereinabove is the wave character of the distribution of the electromagnetic field 
within and around a body. These lead to a highly complicated image reconstruction theory, i.e. 
the problem of diffraction tomography. The linear, ray approach applicable to X-ray 
tomography does not work properly with regard to microwave tomography. Thus, a need exists 
for advanced, non-linear diffraction approaches to the process of image reconstruction. 
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SUMMARY OF THE PRESENT INVENTION 



[0015] The present invention comprises a system and method for non-destructive 
functional imaging and mapping of electrical excitation of biological tissues with the help of 
electromagnetic ("EM") field tomography and spectroscopy using a sensitive material or solution 
introduced into the biological tissue or the circulation system, characterized by having 
dielectrical properties that are a function of electrical field generated by the biological tissue, 
plurality of EM field sources-detectors located around a biological object, so an object imder a 
study is inside an EM field domain, and a control subsystems functionally coupled to the 
plurality of sources-detectors for selectively controlling function of the plurality of sources- 
detectors and for detected EM field from the plurality of sources-detectors so that multiple 
modality EM field is generated on a selected plurality of sources-detectors and detected by a 
selected pluraUty of sources-detectors after interacting with the biological object and a signal 
inversion means operably connected to the control means for inversion of EM fields detected by 
a plurahty of EM field detectors so an image of the biological object and the spread of electrical 
excitation of the biological object are reconstructed. 

[0016] The utility of ttiis invention encompasses many fields of medicine. In 
particular, the embodiments of the present invention relate to non-destructive functional imaging 
and mapping of electrical excitation of biologically excited tissues, including but not limited to 
cardiac tissues, nervous tissue and musculoskeletal tissue. This allows for rapid and accurate 
assessment of functional and pathological conditions of biological tissues and localization of 
sources of irregularities (arrhythmogeneties) in the process of an electrical excitation of 
biological tissues. The identification of the source of such irregularities (arrhythmogeneties) has 
enormous practical utility in selecting a proper treatment and therapy strategy. 

[0017] Further areas of applicability of the present invention will become apparent 
from the detailed description provided hereinafter. It should be understood that the detailed 
description and specific examples, while indicating the preferred embodiment of the invention. 
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are intended for purposes of illustration only and are not intended to limit the scope of the 
invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0018] Further features, embodiments, and advantages of the present invention will 
become apparent from the following detailed description with reference to the drawings, 
wherein: 

[0019] Fig. 1 is a graphical representation of the changes of myocardial s" 

following short time 20%, 40%, 60% and 1 00% blood flow reduction; 

[0020] Figs. 2A & 2B are graphical illustrations of the spectral changes in 

myocardial permittivity e' and resistance p, respectively, during 10% hypoxia; 

[0021] Figs. 3 A & 3B are graphical illustrations of the spectral changes in 
myocardial permittivity s' and resistance p, respectively, during 2 hours acute ischemia; 

[0022] Fig. 4 is a graphical illustration of the changes in myocardial 

dielectric properties (e®) for 2 week old canine myocardial infarction; 

[0023] Figs. 5A & 5B are reconstructed electromagnetic longitudinal 
tomographic images of an excised canine heart together with anatomical slices; 

[0024] Fig. 6 is a schematic view of a method for non-destructive 
functional imaging and mapping of electrical excitation of biological excited tissues 
using electromagnetic tomography and an injection of a sensitive (contrast) material 
(solution) into the biological tissue or in circulation system, characterized by having a 
dielectrical properties that is a function of electrical field, generated by biological excited 
tissue (Case B) compared with traditional method of microwave tomographic imaging 
(Case A); 

[0025] Fig. 7 is a schematic view of an EM field tomographic 

spectroscopic system for non-destructive functional imaging and mapping of electrical 
excitation of biological tissues using a sensitive (contrast) material (solution) injected 
into the biological tissue or in circulation system, characterized by having a dielectrical 
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properties that is a function of electrical field, generated by biological excited tissue, in 
accordance with the preferred embodiments of the present invention; 

[0026] Fig. 8 is a block diagram of one of the EM field source-detector 
clusters of Fig. 7, wherein the cluster is in its source state; 

[0027] Fig. 9 is a block diagram of one of the source-detector modules of 

Fig. 8; 

[0028] Fig. 10 is a block diagram of the R-channel module of Fig. 8; 
[0029] Fig. 1 1 is a block diagram of one of the IF detector clusters of Fig. 

7; 

[0030] Fig. 12 is a block diagram of the control system for the EM field 

clusters and IF detector clusters of Fig. 7; 

[0031] Fig. 13 is a block diagram illustrating the integration of the control 
system of Fig. 1 1 with the system of Fig. 7; 

[0032] Fig. 14 is a block diagram of the EM field soxu-ce-detector cluster 
of Fig. 8, wherein the cluster is in its detector state; 

[0033] Fig. 15A is a flow diagram of the direct problem solver portion of 
the images reconstruction process; 

[0034] Fig. 15B is a flow diagram of the inverse problem solver portion of 
the images reconstruction process; 

[0035] Fig. 15C is a flow diagram of the gradient calculation portion of 

the images reconstruction process; and 

[0036] Fig. 16 is an illustration of an exemplary application of fimctional 
and electrical potential mapping according to a method of the present invention. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0037] Referring now to the drawings, in which like numerals represent like 
components throughout the several views, the preferred embodiments of the present invention 
are next described. The following description of the preferred embodiment(s) is merely 
exemplary in nature and is in no way intended to limit the invention, its application, or uses. 
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[0038] The major idea of the present invention is to use a sensitive material 
(solution) injected into a biological tissue 19 or in the circulation system in order to detect and 
map an electrical excitation of the biological tissue 19. Fig. 6 is a graphical illustration of a 
method for non-destructive functional imaging and mapping of electrical excitation of biological 
excited tissues 19 using electromagnetic tomography and an injection of a sensitive (contrast) 
material (solution) into the biological tissue 19 or in circulation system, characterized by having 
a dielectrical properties that is a function of electrical field, generated by biological excited tissue 
19 (Case B) compared with traditional method of microwave tomographic imaging (Case A). 
The dielectrical properties of this material (solution) are a function of electrical field, generated 
by biological excited tissue 19 itself, so they can be reconstructed with the help of microwave 
(electromagnetic) tomography, and consequentially the spread of electrical excitation of 
biological tissue 19 can be non-invasively reconstructed. 

[0039] Fig. 7 is a schematic view of an EM field tomographic spectroscopic 
system 10 for non-destructive functional imaging and mapping of electrical excitation of 
biological tissues 19 using a sensitive (contrast) material (solution) injected into the biological 
tissue 19 or in circulation system, characterized by having a dielectrical properties that is a 
function of electrical field, generated by biological excited tissue 19, in accordance with the 
preferred embodiments of the present invention. As illustrated in Figs. 7, the system 10 includes 
a working chamber 12, a plurality of "EM field source-detector" clusters 26, an equal number of 
intermediate frequency ("IF") detector clusters 28, and a control system (not shown in Fig. 7, but 
illustrated in block diagram form in Fig. 11). Although only two EM field source-detector 
clusters 26 and two IF detector clusters 28 are shown, it should be clear that a much larger 
number of each, sometimes denoted herein by A^, may (and preferably should) be used. 

[0040] The working chamber 12 may be a watertight vessel of sufficient size to 
accommodate a human body, and may have different shapes and sizes, the selection of which 
would be readily apparent to one of ordinary skill in the art. The working chamber 12 and its 
EM field clusters 26, as well as the IF detector clusters 28, may be mounted on carts in order to 
permit the respective components to be moved if necessary, and the carts may then be locked in 
place to provide stability. 
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[0041] Fig. 8 is a block diagram of one of the N EM field clusters 26 of Fig. 7, 
wherein the cluster 26 is in its source state. Each EM field cluster 26 is a main operation unit 
that may fimction as an electromagnetic field generator (i.e., an electromagnetic source) or as an 
electromagnetic field detector. Each cluster 26 has a plurality of source-detector modules 30, 
one reference channel ("R-channel") module 32 and a pair of distribution blocks 64, 66, as well 
as at least two precision attenuators. The nvunber of source-detector modules 30 (three being 
shown here) in each EM field cluster 26 may sometimes be denoted herein by M In general, the 
more source-detector modules 30 that are used, the greater the precision of the system 10. 
However, because of the large amounts of data that are created, it may be preferable to use 
between 500 and 1500 source-detector modules 30, with an optimum number near 1000, 
organized into N EM field clusters 26, with the value of selected based generally on 
manufacturability and convenience. 

[0042] Fig. 9 is a block diagram of one of the M source-detector modules 30 of 
Fig. 8. Each source-detector 30 includes a BPSK modulator 34, a power amplifier 36, a direct 
uncoupler 38,. a switch 40, a low noise amplifier ("LNA") 42, a mixer 44, a programmable gain 
amplifier ("PGA") 46 and an antenna 48. The switch 40 functions to connect the antenna 48 into 
the system 10 as an EM source or as an EM detector. When connected as a source (i.e., when the 
switch 40 is in the lower of the two positions shown in Fig. 9), an input signal provided by one of 
the distribution blocks 64 (as shown in Fig. 8) is modulated by the BPSK modulator 34, 
ampHfied by the amplifier 36 and uncoupled by the direct uncoupler 38 before passing through 
the switch 40 to the antenna 48. On the other hand, when connected as a detector (i.e., when the 
switch 40 is in the upper of the two positions shown in Fig. 9), the signals received by the 
antenna 48 pass through the switch 40 to the LNA 42 where they are amplified and then mixed 
with a reference signal provided by the second distribution block 66 (as shown in Fig. 8) and 
then ampUfied again by the PGA 46. 

[0043] Fig. 10 is a block diagram of the R-channel module 32 of Fig. 8. As 

described previously, there are preferably a plurality (M) of source-detector modules 30 in each 

EM field cluster 26 but only a single R-channel module 32. The R-channel module 32 includes a 

switch 50, an adder 52, a direct imcoupler 54, an LNA 56, a mixer 58 and a PGA 60. The switch 
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50 controls whether the R-chaiinel module 32 is in its source state or its detector state. When the 
R-channel module 32 is in its source state (i.e., when the switch 40 is in the upper of the two 
positions shown in Fig. 10), output signals from the source-detector modules 30 are passed 
through the adder 52 and the direct uncoupler 54 and are amplified by the LNA 56 before being 
mixed with a reference signal and amplified again by the PGA 60. On the other hand, when the 
R-channel module 32 is in its detector state (i.e., when the switch 40 is in the lower of the two 
positions shown in Fig. 10), a reference signal is passed straight through to the source-detector 
modules 30 where it is coupled with the signals received by the respective antennae 48. 

[0044] Fig. 11 is a block diagram of one of the N IF detector clusters 28 of Fig. 7. 
Each IF detector cluster 28 includes a family of Af+1 digital correlation detectors 70 for Af test 
signals (one from each of the source-detector modules 30 in a corresponding EM field cluster 26) 
and one reference channel signal. These digital detectors 70 allow for the informative/working 
bandwidth of the signal to be selectively passed while restricting other artifacts. Each IF detector 
cluster 28 also includes a cluster manager, a bus, and a power supply. 

[0045] Fig. 12 is a block diagram of the control system for the EM field clusters 
26 and IF detector clusters 28 of Fig. 7. The control system includes a control computer 14, an 
imaging computer 15, a synchronization unit 16, a reference module 18, a distribution network 
20, a calibration appliance 22 and a power supply 24. The control computer 14 controls the 
overall system fimction, data acquisition, system tuning and calibration and transforms all raw 
data to the imaging computer 15 for fiuther data inversion and imaging. The control computer 
14 may be a conventional personal computer, such as an Intel-based advanced-version PC, with 
an RS-488.2 port and appropriate software to control the system 10. The synchronization unit 16 
is a module that includes a system manager and a system hub. Together, they provide data 
exchange with the control computer 14 (preferably via a USB 2.0 or Firewire link) and the 
control managers of the various clusters 26, 28, and also provide synchronization of system 
operations. 

[0046] The reference module 18 includes two generators, one or more thermostats 

for temperature stabilization of the fiinction of the reference channels, a BPSK modulator for 
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phase-modulation, power dividers, attenuators and the Hke. The two generators are precision 
generators that generate stable CW signals: Carrierref and LOref. These generators are controlled 
and tuned by the control computer 14 through an interface. The distribution network 20 is a 
commutation unit for receiving the carrier and local oscillator reference signals (Carrierref and 
LOref) and the Rr and Rtr reference signals (Rxref and Rtrref) from the reference module 18 and 
distributing them to each of the EM field clusters 28. 

[0047] The calibration appliance 22 is used for calibration and fine-tuning of the 

system 10. The calibration appliance 22 includes a calibration source, one or more (preferably 
two) calibration antennae, precision drives and one or more (preferably three) calibrated 
phantoms. Calibration antennae and phantoms may be precisely positioned at any point inside 
the working chamber with the help of precision positioning drivers. The isolated power supply 
24 provides stable power for the system. One power supply suitable for use with the present 
invention is a 190/380 3 -phase, 10 kVA AC network power supply. Of course, the exact 
requirements for the power supply 24 may depend upon the power system specifications of the 
country in which the system 10 is to be operated. 

[0048] Fig. 13 is a block diagram illustrating the integration of the control system 
of Fig. 11 with the system 10 of Fig. 7. Each EM field cluster 26 is disposed adjacent the 
working chamber 12 such that its antennae 48 are located on or near the surface of the chamber 
12. The outputs of the source-detector modules 30 and the R-channel module 32 of each EM 
field cluster 26 are connected to a corresponding IF detector cluster 28, and each IF detector 
cluster 28 is connected to both the corresponding EM field cluster 26 and the synchronization 
unit 16. The inputs of each EM field cluster 26 are connected to the distribution network 20. 
The distribution network 20 includes at least four distribution blocks 68, which may be 34- 
channel power dividers, and a system bus for distributing the various reference signals (Carrierref, 
LOref, Rfref and RtTref) to the EM field clusters 26. As illustrated in Fig. 12, one set of the four 
signals is provided to each EM field cluster 26. These signals are denoted Carrierj, LOi, Rrj and 
Rtri, where the first EM field cluster 26 (i.e., the one shown in Fig. 2F) receives Carrier], LOi, 
Rri and Rtri, the second EM field cluster 26 (not shown) receives Carrier2, LO2, Rr2 and Rtr2, 
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and so forth. Finally, as described previously, Carrierref and LOref are provided to the distribution 
network 20 by the reference module 18. ' 

[0049] In use, the working chamber 12 is filled with one of a variety of solutions 
or gels 17 selected to match and provide biological compatibility with a biological tissue object 
19 to be studied. Suitable solutions 17 may include, but are not limited to, water, salt solutions, 
sugar solutions, fatty emulsions and the like; these solutions may also be used as gel 
components. The object 19 to be studied may be injected with a sensitive material (solution) (or 
distributed in the object 19 via the circulation system) whose dielectrical properties are a 
function of the electrical field generated by the biological excited tissue 19 itself, so that they can 
be reconstructed via microwave (electromagnetic) tomography. Preferably, the injection 
materials or solutions are a multicomponent media that includes ferroelectrics, such as barium 
modified strontium titanium oxide, of different grain sizes ranging from 0.5-100 |xm. The grains 
preferably also include different shapes, including spheres, ellipsoids and cylinders. The 
materials may also include selected potentiometric dyes, such as merocyanine, rhodamine, 
cyanine, oxonol and naphthyl styryl, and/or selected potentiometric liquid crystals, such as 
MBBA, 7CB. 

[0050] After injection, the object 19 is then positioned inside the working 
chamber 12 and the system 10 is activated. During operation of the system 10, each Carrieri 
signal fi-om the signal generator in the reference module 18 is provided to a source-detector 
module 30, operating in its source mode as shown in Figs. 8 and 9, where it is modulated using 
phase-shift modulation (in case of phase characterization) by pseudo-random code in order to 
distinguish each transmitting antenna 48 or source fi"om the other antennae/sources 48, which are 
transmitting simultaneously. As described previously, the resultant signal is next amplified 
before passing through the direct uncoupler 38 to the appropriate source antenna 48. As a result, 
an incident EM field ("Ejnc")? corresponding to the respective antenna 48 or channel, is formed in 
the vicinity of the object 19 under study. In addition, part of the signal creating the Einc field is 
uncoupled and passes to a receiver in the R-channel module 32 (one for each EM field cluster 
26). In the R-channel module 32 this signal is mixed with a reference signal Rri. By 
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subsequently comparing the resultant output with a known signal, Ejnc may thus be determined 
precisely as described below. 

[0051] After interacting with the object 19 of interest, each "interferenced" or 
scattered EM field (Esct) is detected by a corresponding detecting antenna 48 operating in its 
detector mode. Fig, 14 is a block diagram of the EM field source-detector cluster 26 of Fig. 8, 
wherein the cluster is in its detector state. The same reference signal Rrj described in the 
preceding paragraph is injected into the source-detectors 30 of the EM field cluster. 26 (operating 
in detector mode) immediately downstream from the detecting (receiving) antenna 48. This 
allows for the R-channel signal Rrj, which is known precisely, to pass through all parts of the 
detector 30 through which the Esct signal is passed. Therefore, an injection of the R-channel 
signal into the measuring portions of the source-detectors 30 in both source and detection mode 
allows for a significant decrease in artifacts caused by temperature and temporary instabihty of 
the channel electronics. 

[0052] The data and other information gathered by the system 10 is provided to 
the imaging computer 15. A novel process is carried out by the imaging computer 15 to solve an 
inverse problem of electromagnetic field tomography. It is based on the solution of a non- 
simplified three-dimensional ("3D") vector using Maxwell's equations. Figs. 15 A, 15B and 15C 
are flowcharts of this process. The method uses an iterative procedure based on gradient 
calculation approach with the following features, among others: 

• The method is based on minimization of the difference between model scattered fields 
and measured scattered fields. 

• The method uses the Tichonov's type of regularization. 

• One type of the calculation mesh is used in the method. 

• One step of the iterative procedure is performed as solving of the two sets of direct 
problems of the same dimension: modehng of the so-called direct wave and modeling of the 
inverse wave. 

• Both the direct wave and the inverse wave are calculated using nonreflecting boundary 
conditions. 

• Both the direct wave and the inverse wave are calculated on the same rectangular mesh. 
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• In order to solve the direct problem a conjugate gradient method ("CGM") is used. 

• One step of the CGM uses the sine Fourier transform. 

• The wave equation for non-uniform media is used to solve the direct problem. 

[0053] From a mathematical point of view, the methodology utilized in EM field 
tomography is an inverse problem. It may be formulated in terms of complex dielectrical 
properties s and electrical and magnetic fields -E,H , The basis is a set of the Maxwell's 
equations: 



. 4;r . i€0) _ 
c c 

cirlE = -—H (1) 
c 

div{sE) = 0 
divH = 0 



where represent electrical and magnetic fields, respectively, and all other notations are 
standard. 

[0054] It is more practical to rewrite these equations in a form of non-uniform 
wave equations: 



[0055] 



V'E + k'E- divgradE = 0 (2) 



[0056] 



where = 



2^ 



s , and X is a wavelength in vacuum. The EM field 



tomographic system could be schematically represented as a chamber with the set of antennae on 
the surface of the chamber. As described previously, some antennae function as EM field 
sources while the others fimction as EM field detectors. It is usefiil to divide electric field E 



into incident E^ field and scattered fields, : 



(3) 
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TTiP pniiation (2) can be rewritten in 
„here,- Utoenumbar of a particular «iner or sour«. KeequanonC ) 

the form: 

V'E, +k%j -graddivE, =(kl -e)E,j (4) 
„.e. , I a wave for and . « pro^oM ^ 

antenna number J . 

,0057, A.obiect™aybedescH..dasadis.ribu,ionotaie.ec«ica:pe— in 

the volume. 

(0058, A receiver antenna records ^ signa,, which reflec. boU. incident and 
scattered fields. 

(A\ wp need to use some boundary 
mni,Ql In order to solve the equation (4) we need to 
[0059] in oraer . ^ • we are using nonreflecting boundary 

solved as a separate problem. 

[0060, m ...^ u,i„„ tt,e FDTD method we have found that 

edge of the antenna: 

£(,)=ff[(n'xE(r'))xV'GH-(«'-£)V'G]d^' (5) 

where the Green function is: 



E,(r') = EoCOs{ny/a) 
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The equation (5) shows good agreement with the experimental results. The same type of antenna 
may be used as a receiver. To describe the process of a signal recording we use the reciprocity 
principle. It could provide different expressions for recorded signal. We use two of them in our 
calculations. First: 

Sy=qlllE,{r)Ej(r)A^^^^^ (6) 

where S^j is the signal received by antenna number lin the situation where antenna number 

7 works as a transmitter or source, C, is a constant, E. is an electrical field distribution 

produced by the detecting antenna, Ej is an electrical field distribution produced by the source 

antenna, Ae is dielectrical permittivity distribution in the object, and the integral is taken in the 
domain V where the object is located. The equation (6) is used in the inverse problem solution. 

S,=C,j^[E,xHj-EjxH,]ds (7) 

5 

where the integral is taken over the surface of the domain. 

[0061] Direct problem solver . In order to solve direct problems we use a 
conjugate gradient method with a preconditioner. In order to do that, equation (4) should be 
rewritten in the form: 

AE, =V'E^^klE^-graddivE^ =i^kl^k')E^ ^(^ -k')E, (8) 

where k^^ is an average value of k . The preconditioner operator can be constructed as a first 
step of the iterative process: 

^r' = (kl -k')E: +(k', -e)E, (9) 

Taking into account the fact that the left side of equation (8) is an expression with constant 
coefficient, (9) can be solved at step 1575 using sine-type Fourier transform for the case with 
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zero boundary conditions on the bound of calculation domain. Then R.A. James's method 
(originally invented for static problems, but subsequently developed* for electromagnetic 
problems) is applied to make boundary conditions nonreflected. This technique creates a very 
robust and effective method. Computational experiments show that the iterative process appears 
to work with any reasonable contrasts and provides nonreflecting conditions with very high 
accuracy. Using a sine-type Fourier transform at step 1575 can make calculations 8 times faster 
than with the regular Fourier approach. 

[0062] Fig. 15A is a flow diagram of the direct problem solver portion 1535 of 
the image reconstruction process. The direct solver 1535 is used only for inverse problem 
solving. The input data in this case is the dielectric properties distribution in the form of a 3D 
array, which is received at step 1560. For the first step of the iteration, this input data is received 
from external input, while in subsequent iterations it is received from the previous iteration. 
Next to occur, at step 1565, is the preparation of the parameters and arrays, which do not change 
during the direct problem solving process: the wave number, the computational grids, and the 
Green function for the uniform space. After that, the iterative procedure of the conjugate 
gradients takes place at steps 1570-1580. First, the source member of equation (4) is calculated 
at step 1570. Then, every step of the conjugate gradient method requires fast Fourier transforms 
of the source functions, as shown at step 1575. hi order to stop iterations the convergence of the 
process is checked at step 1580. Once the iterative procedure is finished, the non-reflecting 
boundary conditions have to be implemented at step 1585. Finally, the output of the process 
1535 is created at step 1590. The output comprises arrays containing the electric fields inside of 
the computational domain and signals on the receivers for all transmitter positions. 

[0063] Inverse problem solver . In order to solve the inverse problem in 
microwave tomography we apply the gradient method. In the case of a three-dimensional vector 
in cylindrical geometry this method needs significant modifications when compared with two 
dimensional and scalar cases. In general the inverse problem in EM field tomography can be 
formulated as a minimization problem: 
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Ae] = ^\Sy"" -S;'"" f +an[s] (10) 

where Sy^""' are the theoretical values of the signal, Sy^^^' are experimental values of the signal, 
and the last term is the Tichonov's regularization functional, 

[0064] The key point of any minimization procedure is the method of a gradient 
calculation. It was proven that the gradient of functional in our case is: 

y'[f] = X(^yG,/(5f- -5;^-)+an' (ii) 

ij 

where Ej and Gy are solutions of the following equations: 

V^Ej + k^Ej - V{VEj ) = Fj (12) 
V%+k%-S/(yG,) = P, (13) 

Functions F, and describe the field pattems for antennae 48 being used as sources and 
detectors, respectively. 

[0065] Direct computation using the equation (11) is very time consuming even 
in the 2D case and cannot be effectively applied in the 3D case. The reason is that every step 
requires NxM number of direct problems to be solved, where N is the number of transmitters, 
and M is the number of receivers. It was shown in our previous "scalar" work, and can be 
generalized in the vector case, that the function 

Zj=TGy(st"-sr'y (14) 

can be the solution of the following equation: 

S7'Zj+k'Zj - V(VZ,) = ^P,(S^^' -S^'') (15) 
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This makes it necessary to solve only two direct problems on each iterative step. 

[0066] The calculation of the sum in the right side of equation (15) continues to 

be a difficult problem, because it requires summation on all receivers for all cells of the 
computational mesh. In order to overcome this obstacle, a two-step procedure may be applied. 
First, the following may be calculated on the surface of the computational domain: 

^;=z p,(^t""-sr') (16) 

This needs significantly less computational effort compared to the calculation of the right part of 
equation (1 5). Second, the following equation may be solved with those boundary conditions: 

V'Z^'+^o'Z^' =0 (17) 

Equation (17) is the equation with constant coefficients and can be easily solved using sine-type 
FFT. 

[0067] Finally, one step of the gradient method procedure requires solving two 

direct problems (equations (12) and (15)) plus one equation (equation (17)) with constant 
coefficients. In general this procedure looks to be the fastest known in literature. 

[0068] One step of the iterative procedure can be implemented as: 

=^"-J'[8"]5 (18) 

where an iterative step is chosen in a trial method. The limitations on the upper and lower 
bounds of the values of the dielectric properties and the values of the dielectric properties on the 
bound of the object are applied in this step. 

[0069] Fig. 15B is a flow diagram of the inverse problem solver portion 1500 of 
the image reconstruction process. At step 1505, the input data is received. The input data for the 
inverse problem solver 1500 includes physical and geometrical parameters of the computational 
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process: the sixes of the computational domain, the working frequency, the maximum number of 
iterations and the signals from the antennae 48. Next to occur, at step 1510, is the preparation of 
the parameters and arrays, which do not change during the inverse problem solving iteration 
process: the wave number, the computational grids, and the Green function for the uniform 
space. After that, the iterative procedure of calculating the gradient of the residual function 
(equation (11)) itself takes place at steps 1515-1520, including the gradient calculation process 
itself at step 1515, In order to stop iterations the convergence of the process is checked at step 
1520. This involves comparing the value of the residual error with the estimated experimental 
error. Once the iterative procedure is finished, the non-reflecting boundary conditions have to be 
implemented at step 1525. Finally, the output of the process 1500 is created at step 1530. The 
output comprises the dielectric properties distribution in the form of a 3D array. 

[0070] Fig. 15C is a flow diagram of the gradient calculation portion 1515 of the 

image reconstruction process. The direct wave is calculated at step 1535 according to equation 
(12), followed at step 1540 by the calculation according to equation (16) of the source for back- 
propagating wave on the bounds of the computational domain. Then, at step 1545, the source of 
the back-propagating wave is calculated in the volume of the computational domain according to 
equation (17), and the back-propagating wave is calculated by solving equation (13) at step 1550. 
Finally, the gradient is calculated according to equation (1 1) at step 1555. 

[0071] The image reconstruction algorithm of this invention includes a number of 

benefits. For example, using the nonreflecting boundary conditions plus sine-type FFT makes 
the direct problem solver of the invention the most effective one. Further, the proposed way to 
calculate the so-called back wave (equations (15), (16), (17)) allows working in real 3D multi 
antenna configuration. In addition, the method of signal calculation (equation (7)) is 
distinguished from any others and allows simulating the work of each antenna with high 
precision, and the mathematical algorithm itself is essentially parallel, which is particularly 
advantageous for parallel computing. 

[0072] Fig. 16 is an illustration of an exemplary application of functional and 

electrical potential mapping according to a method of the present invention. As disclosed above, 
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the invented system and method allows for both functional imaging of biological objects 19 and 
mapping of electrical excitation of biological tissues 19. This unique feature of the invention can 
be achieved as a two-step imaging process, as illustrated in Fig. 16. First, the system 10 is used 
as an electromagnetic tomographic spectroscopy tool for non-invasive assessment of functional 
and pathological conditions of a biological object 19 and the location of any areas of diseased 
tissues 19. For example, in Fig. 16, cardiac tissue with suspicious areas of myocardium is 
functionally imaged in Step 1 to identify areas of ischemia/infarction. Second, if the location of 
a diseased target is within a biologically excited tissue 19 (such as cardiac tissue, nervous tissue, 
musculoskeletal tissue, or the like), then a dielectrical contrast (sensitive) substance (solution) 
may be injected into the tissue (circulation system), thus allowing for non-invasive mapping of 
electrical excitation and the location of a source of irregularity (arrhythmogenety). Since the 
functional/pathological condition of the tissue in and near this source is known from the first step 
of the EM field imaging, it allows for imique and extremely valuable information to be 
developed for use in determining tissue viability and in choosing a strategy of further treatment 
and therapy. For example, in Fig. 16, electrical potentials in the areas of ischemia/infarction in 
the cardiac tissue are mapped to identify possible sources of arrhythmia and tissue functional 
viability. Further, if an ablation of this tissue is chosen than the same system can be used and 
electromagnetic energy can be focused precisely into the target area according to known 
techniques such as those described in the aforementioned U.S. Patent Nos. 5,715,819, 6,026,173 
and 6,333,087. 

[0073] Of course, it should be apparent that it may not be necessary for the first 
step to be conducted at all. For example, other, more conventional, tests or diagnostic tools may 
alternatively be utilized to identify a particular tissue or area of tissue that are believed to warrant 
closer study. The second step described above may obviously be used to map the electrical 
excitation of the tissue and to isolate irregularities in the manner described herein. 

[0074] The invention uses multimodality EM field: multi-frequency as well as 

multi-polarization with special timely fashion, synchronized with a biological electrical signal. 

This allows for both a more efficient and accurate image reconstruction algorithm and a 

combination of imaging and tissue spectroscopy for better assessment of tissue viability. The 
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code-division technology may be used to distinguish each particular EM field source antennae 
from a plurality of simultaneously operating antennae. The multi-frequency EM field is 
preferably within a frequency range of about 50 KHz to 10 GHz and may be organized in a 
number of clusters {M) with closer frequencies in each frequency cluster. In this case, using N 
close frequencies (with differences in about lOMHz) allows for more effective images 
reconstruction having ahnost frequency independent tissue dielectrical properties. Since the 
reconstruction algorithm described above deals with so-called scattered EM field it requires 
measurements of the distribution of EM field in the working chamber 12 filled with a matching 
solution 17 but without an object 19 - i.e., a so called "empty" field. This multi-frequency pack 
of radiation (Aldose frequencies) may be used instead of additional measurements of an "empty" 
field within the working chamber 12. This helps to solve stabihty problem of the system 10 and, 
finally, to improve reconstructed images. The M frequency clusters varied significantly (about 
or more than 50-100 MHz) are needed for tissue spectroscopy purposes. 

[0075] In a fiirther feature of the present invention, the fiinctional viability of 
blood vessels may be assessed by introducing into the circulation system of the biological tissue, 
through injection or otherwise, a dielectrical contrast solution that is characterized by having 
dielectrical properties significantly different from those of the blood normally carried by the 
blood vessels. The system and methods described above may then be applied. Using code- 
division technology the total data acquisition time of the system is decreased to 5-10 msec. 
Clinically-proven iodine-based radiopaque agents, including but not limited to diatrzoate 
meglumine, or intralipid or other solutions and their mixtures may be used for this purpose. 

[0076] Based on the foregoing information, it is readily understood by those 

persons skilled in the art that the present invention is susceptible of broad utility and application. 

Many embodiments and adaptations of the present invention other than those specifically 

described herein, as well as many variations, modifications, and equivalent arrangements, will be 

apparent from or reasonably suggested by the present invention and the foregoing descriptions 

thereof, without departing from the substance or scope of the present invention. Accordingly, 

while the present invention has been described herein in detail in relation to its preferred 

embodiment, it is to be imderstood that this disclosure is only illustrative and exemplary of the 
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present invention and is made merely for the purpose of providing a full and enabling disclosure 
of the invention. The foregoing disclosure is not intended to be construed to limit the present 
invention or otherwise exclude any such other embodiments, adaptations, variations, 
modifications or equivalent arrangements; the present invention being limited only by the claims 
appended hereto and the equivalents thereof Although specific terms are employed herein, they 
are used in a generic and descriptive sense only and not for the purpose of limitation. 
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